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Abstract 

We investigate the two-photon corrections to the process pp — t- e'^e" at large momentum trans- 
fer, aimed to access the time-like nucleon form factors. We estimate the two-photon corrections 
using a hard rescattering mechanism, which has aheady been used to calculate the corresponding 
corrections to elastic electron-proton scattering. Using different nucleon distribution amplitudes, 
we find that the two-photon corrections to the pp — )• e"'"e~ cross sections in the momentum transfer 
range 5-30 GeV^ is below the 1 % level. 

PACS numbers: 
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I. INTRODUCTION 



Electromagnetic form factors (FFs) provide important information on the structure of 
the nucleon. Consequently, there has been much effort in their measurement both in the 
space-like as well as in the time-like regions. 

The space-like electromagnetic FFs, which provide information on spatial distributions 
of quarks in the nucleon, can be investigated in elastic electron-proton scattering; for recent 
reviews see e.g. Refs. jl-^. Two experimental methods exist to extract the ratio of electric 
(Ge) to magnetic [G^j) proton FFs. Historically, the first method involves unpolarized 
measurements employing the Rosenbluth separation technique, which gives a direct acess to 
the space-like FFs through the slope and intercept of the e-dependence of the cross section 
in the one-photon (I7) exchange approximation : 

rfa = C(g^e)te(Q^) + -G|(Q2)l, (1) 
L r J 

where £ and are the virtual photon polarization parameter and virtuality respectively, 
r = Q'^/Am%, with m^v the nucleon mass, and where C is a known phase space factor. In 
more recent years, polarization experiments using polarized electron beams on polarized 
targets or measuring the recoil nucleon polarization, in the elastic e—p scattering, provided 
another experimental way to access Ge/Gm- The ratio of polarization of the recoiling proton 
perpendicular to its motion (Pj) to polarization along its motion (Pi) is directly related to 
the ratio of electric to magnetic proton FFs 

Pt I 2e GEm 



Pi y T{l + e) GMiQ^y 



(2) 



Such polarization experiments have been performed for momentum transfers up to 8.5 GeV^ 
to date and have shown that the ratio of the electric to magnetic proton FFs is decreasing 
with increasing momentum transfer j4-7|. This finding is in contrast to the well known 
scaling-behavior of Ge/Gm determined by the Rosenbluth separation technique. The dis- 
crepancy between data of unpolarized Rosenbluth measurements and polarization experi- 
ments has t rigger ed a whole new field studying the influence of two-photon (27) exchange 
corrections 8-13|; see Ref. [l^ for a recent review and references therein. The finding of 
those works is that two-photon exchange corrections to the Rosenbluth cross section are a 
possible explanation for the discrepancy, whereas the 27 exchange effects do not impact the 



polarization transfer extraction of Ge/Gm in a significant way. Recently, a first empirical 
extraction of the three 27-exchange amplitudes to elastic electron-proton scattering has been 
performed [l^ based on measurements of cross sections [l^ and polarization observables 
at a common value of four-momentum transfer, around = 2.5 GeV^. It confirms that 
a common description of unpolarized measurements and of polarization observables invokes 
empirical 27-amplitudes with relative magnitude up to about 3 %. 

The measurements of nucleon FFs at space-like momentum transfers, through elastic 
electron-nucleon scattering, are complemented by measurements in the time-like region, 
through the crossed processes pp — ?■ e~^e~ and e~^e~ — NN, which access the vector mesonic 
excitation spectrum of hadrons. The latter process has been measured in recent years at 
e~^e~ facilities, such as Da$ne, CLEO, and BABAR. These measurements have revealed that 
the nucleon FFs at time-like momentum transfers are significantly larger than their space- 
like counterparts when considering the same magnitude for the virtuality = —Q^. In 
particular, for momentum transfers with magnitude around 10 GeV^, the time-like FFs were 
found to be enhanced by a factor of two. New measurements are planned in the near future at 
BES-III and at PANDA@FAIR, bringing time-like (positive) values around 20 GeV^ into 
reach. Such new measurements will explore the at present still largely uncharted time-like 
region in much greater detail and complement our picture of the nucleon. 

The time-like FFs are complex quantities due to the interactions of the hadrons in the ini- 
tial and final state, respectively. Their absolute values can be determined from measurements 
of the angular distribution of the unpolarized cm. cross section in the 17-approximation : 



|Gm|2(1 + cos^^) + -\GE\^sm^9 

T 



(3) 



whereas the phases are related to polarization observables. Since 27-exchange plays a cru- 
cial role in the extraction of electromagnetic FFs in the space-like region, investigating its 
influence in the time-like region seems to be an obvious task. Even though some theoretical 



works have been done 



18 



19| . there are no comparable calculations so far to estimate the 



27-exchange corrections for the corresponding time-like processes. 

In this work we investigate the 27-exchange corrections to the process pp — e^e^ at 
large momentum transfer g^. To provide a first estimate of the corrections we consider a 
perturbative QCD (pQCD) factorization approach, which has already been used to calculate 
the corresponding corrections to elastic electron-proton scattering jl3| . 
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The paper is organized as follows: The general formalism including 27-exchange is pre- 
sented in Section [Tll In Section IIIII we estimate the hard 27-exchange contribution at large 
momentum transfers by relating the 27-exchange amplitude to the leading twist nucleon 
distributions. The results of the calculation are discussed in Section IIVI Some concluding 
remarks are given in Section |Vl 



II. GENERAL EXPRESSION OF THE OBSERVABLES INLCUDING 27- 
EXCHANGE 

In order to describe the annihilation of a proton and an antiproton into a lepton pair, 

p{Pi, ^Ni) + P{P2, AatJ r{ki, hi) + t{k2, h2), (4) 

where X^i, Aatj, hi and /i2 are helicities of the nucleons and leptons respectively, we adopt 
the definitions 

P = — 2 — ' " — 2 — ' ^ = (^'i + Pv . (5) 

and the Mandelstam variables 

s = g' = t = {pi-k2)^ u = {j)i-kif. (6) 

The process can be described by two independent kinematical invariants, which we choose 
as the variables and t. 

The amplitude of the reaction is related by crossing to the corresponding scattering 
amplitude for elastic electron-proton scattering. Neglecting the lepton masses, the matrix 
element including multi-photon exchange is parameterized by three generalized form factors. 
Several equivalent representations exist. Here we use the representation, which was first 



introduced in Ref. 



T = u{k2, h)'y^v{ki, -h) x v{p2, Xn2 



. The matrix element can be written in the form 

u{pi,Xni] 



GmY - F2 —P'' + Fs \P^$ 
niN my 



(7) 

where Gm, F2 and F3 are complex functions of and t. Neglecting the lepton masses 
implicates that the outgoing electron and positron have opposite helicities. 
In the following, we also use the generalized form factor 



In the Born approximation Gm and Ge reduce to the usual proton FFs and do not depend 
on t, while F3 vanishes. In order to identify the I7 and 27-exchange contributions, we 
introduce the decompositions 

GM{q\t) = GM{q') + SGMiq\t), 
GE{q\t) = GE{q^) + 5GE{q\t), 

F,{q\t) ^ 6F,{q\t). (9) 

Gm and Ge are the time-like proton magnetic and electric FFs and -F3, 5Gm and 5Ge are 
amplitudes of order e^, which originate from processes involving the exchange of at least two 
photons. 

To compute the differential cross section of the reaction, we use the center-of-mass (cm.) 
frame, where the momenta of the incoming proton and antiproton have opposite directions. 
In this frame the variable t can be related to the c.m.-scattering angle 9 between the incident 
proton and the outgoing electron. Calculating the cross section up to next order of leads 
to the expression 



\Gm\'^{ 



cos 



+ 2Re[GM5GA/*](l + cos2 



-\GE?smH 



2-Re[GE5GE ] sin^ 9 



2{Re[GMF-i] - ^Re[G£F3*]) Vr(r - 1) cos ^sin^^ 



with 



r 



q' 



C{q' 



T 



(10) 



(11) 



In the 17-exchange approximation, only the first two terms of Eq. ( ITOll contribute to the 
cross section and it reduces to the well known formula of the unpolarized cross section: 

1 



da, 



cm. ,17 



C{q' 



\Gm\\1 + cos' 



IGRMsin^e 



(12) 



The other part of Eq. ( ITOI) represents the interference of I7 and 27-exchange processes. 

In order to determine the imaginary part of the time-like form factors it is necessary to 
study polarization observables. An observable which gives acess to the imaginary part of 
the electric and magnetic form factor is the single spin asymmetry when either the proton or 
antiproton is polarized normal to the scattering plane, which does not require polarization of 



the leptons in the final state. Polarization of the proton or antiproton along or perpendicular 
to its motion, but in the scattering plane, in contrast also requires a polarized lepton. 
The single spin asymmetry can be defined as 

_ da^ - da^ 

^'~da^ + da^' ^^^^ 
where da"^ {da^) denotes the cross section for an incoming nucleon with positiv (negativ) 
perpendicular polarization. In the case of a polarized proton the single spin asymmetry up 
to next order in can be obtained as 

Ay = ^2sin^i (lm[GEGlj]+lm[GE5G*M]+lra[5GEGlj])cose 

dac.m. ^ 

+ Vr(r - 1) {\-m\GMFl\ cos^ B + \^\GeFI\ sin^ | . (14) 

In contrast to space-like processes the single spin asymmetry Ay in the time-like region does 
not vanish in the Born approximation. 



III. CALCULATION OF THE 27-EXCHANGE CONTRIBUTION AT LARGE 

To calculate the 27-exchange corrections in — )■ e^e~ at large momentum transfers 
we consider a factorization approach using the concept of hadron distribution amplitudes 
(DAs). We follow the experience gained by the space- like process ep — )■ ep, for which the 



amplitudes bGu and F3 were computed 
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13| at large momentum transfer Q in the form 



of a convolution of a hard kernel if, which can be calculated in QCD perturbation theory, 
and the nonperturbative contributions which can be related to the DAs of proton and 
antiproton, for instance: 

bGu{(f.e)-^^H[Q\e)^^, (15) 

where the asterisk denotes the convolutions with respect to the participating quark mo- 
mentum fractions. This result represents the leading order contribution with respect to an 
expansion in l/Q. The important feature of such an approach is that the virtualities of both 
photons must be large: —q\ ~ — g| ~ . The corresponding hard subprocess involves only 
one hard gluon exchange and is therefore suppressed by the strong coupling ~ a^. As all 
spectator quarks are involved in the hard scattering process described by Eq. (fT5|) . we shall 
refer to it as the hard rescattering contribution. 
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This mechanism can be simply generahzed to the crossing channel pp — )■ e^e~ at large 
momentum transfer > 0. A typical diagram of the leading pQCD contribution to the 27- 
exchange correction to the annihilation amplitude is illustrated in FigJTJ The simple analysis 



allows to conclude that a description of the corresponding hard time-like subprocess can be 
obtained directly from the space-like one using the crossing symmetry. Therefore we obtain 
that the leading order asymptotic behavior of the time-like 27-exchange amplitudes can be 
represented by the same form as Eq. f|T5|) for the space-like process. 

In the case of g^,t ^ the momenta of proton and antiproton in the c.m.-frame can 
be expressed by two light-like vectors n, n: 




FIG. 1: Diagram for — t- e"'"e including the exchange of two hard photons 




(16) 



The lepton momenta are defined as 




(17) 



where, at large momentum transfer, 77, rj and k±_ can be determined from 



with the obvious restriction < 77 < 1. Note that the kinematic variable 77 can be expressed 
in terms the e~ cm. angle 9 as : 

r/^i(l + cos^) (19) 
The proton matrix element at leading twist level is described by twist-three nucleon DAs 

as: 

4 (0 \e'^''u^W'[\in]upW^[\2n]d„W^[\^n] \pi) = j Dxi 6-'^'+^^''^''^^^ ^p„{xi), (20) 
with measure given by Dxi = dxidx2dx36{l — xi — X2 — X3), and where 

qaW[x] = qa{x)P exp !^ig j dt {n ■ A) {x + tn)^ . (21) 
23| . the function "^aiSaixi) can be expressed as : 

+ T(x.) pi+ [i#7± [7^75iVl , , (22) 

where = represents the large component of the nucleon spinor, C is charge con- 
jugation matrix: C~^7^C = — 7J, and the scalar functions A, V, T stand for the nucleon 
DAs. 

The hard rescattering contribution to the time-like 27-exchange amplitudes SGm, and 



Following Ref. 



s/m\F^ can be obtained from the results for elastic ep-scattering [13| using crossing relations 
for the hard perturbative subprocess. In our kinematics, expressed by Eqs. ([T6 | [T7 t [T8|l . this 
leads to the following substitution: 

Q2 ^ _g2 -^e, V, (23) 



where ( is the kinematical parameter introduced in Ref. 13|. We then obtain for the time-like 
27-exchange amplitudes : 

6GMiq\v) = ( - ) / mv-i)^2yMyuX.) 

J yi2/22/2 X1X2X2 [x2?7 + y2'n - X2y2\ [x2'n + y2'n - ^2y2\ 





(2-K 















F^{q^ 7]) = '"'^"^'^^ l — \ I ^[^^] ^[^»] 2{x2y2 + X2y2)^{yi,Xi) ^^5^ 



q^ \6 J J yiy2y2 X1X2X2 [X2?7 + y2'n - X2y2\ [X2t1 + ^2?? - X2y2 
where $ denotes the specific combination of the nucleon distribution amplitudes: 

^{y^.Xi) = Qu^ [{V + A){V + A) + AT'T] (3,2,1) 

+ QuQd [{V + A')iy + A)+ AT'T] (1, 2, 3) + 2Q„g, [V'V + A' A] (1, 3, 2), (26) 



8 



and the numbers in the brackets define the order of the momentum fractions in the arguments 
of the DAs: V'V{3,2, 1) = V {y^, y2, yi)V {x^, X2, xi) . We also introduced the quark charges 
Qu = +2/3, Qd = —1/3, the fine structure couphng = e^/(47r), and the QCD couphng 
as- 

In general, the time-like amplitudes are complex functions. At tree level, the expressions 
of Eqs. ( 124|) and (!25|) do not contain an imaginary part explicitly. This can be simply 
understood : the s-channel cut requires the on-shell photons (see e.g. diagram in Fig{T]) 
but at large their propagators are highly virtual and hence the tree amplitudes are real. 
Therefore we can obtain nontrivial imaginary contributions only from the loop corrections. 
In particular, computing leading logarithms associated with the renormalization of DAs and 
QCD coupling one obtains imaginary contributions generated by time-like logarithms: 
ln[— — ie] = ln[g^] — in. Such effects can be easily accounted for by using the well known 
formula for the analytic continuation of {20 1: 

«.(-^^) = ^-f^4^ + -, (27) 

1 - z/3oas(r)/4 

where Pq = 11 — 2/3ny is the first term of the /3-function. Eq. (1271) includes resummed 
large corrections ~ PqOs which can be important at intermediate energies where is not 
too small. Similarly, solving the renormalization group equation, one obtains an imaginary 
part originating from the evolution of DAs. However, the resulting imaginary contributions 
provide quite small numerical effects for the regions of which we are going to discuss 
below, see e.g. Ref. [21|. 

As can be seen from Eqs. flMf [25]) the leading behavior of the amplitudes SGm and 
s/m'lfFs goes as 1/g^, whereas 6F2 is suppressed in the large momentum transfer limit, 
since it behaves as We may expect that at intermediate energies ~ 5 — 10 GeV^ the 

effective scale defining the applicability of the perturbative expansion is already large enough 
in order to apply the present formalism. In what follows, we assume that the scale of the 
strong coupling in Eqs. fl2^ [25]) is of order ^ 0.6 g^. 

To evaluate the convolution integrals give n in Eqs. (!24t [251) . we need to consider a model 
description for the twist-3 DAs. In Ref. 22] the asymptotic behavior of the DAs and their 
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In 

(10-3 GeV2) 


r_ 


r+ 


COZ \2£ 


5.0 ±0.5 


4.0 ± 1.5 


1.1 ±0.3 


BLW [25] 


5.0 ±0.5 


1.37 


0.35 


QCDSF [26] 


3.23 


1.06 


0.33 



TABLE I: Parameters entering the proton DA (at /i = 1 GeV) for three parameterizations (COZ, 
BLW, and the lattice evaluation from QCDSF) used in this work. 

first conformal moments can be found as 



A{xi) 
T(x,) 



120XiX2X3/Ar [1 ± r+(l - 8X3)] 

120x1X2X3 fN r_(x2 - Xi), 



120X1X2X3/ 



N 



1 ± - (r_ 
2 ^ 



3x3) 



(28) 



where the DAs depend on the three parameters Jn, r+ and r_. For our calculations, we 
consider two phenomenolo gica l models for the DAs, which have been discussed in the litera- 



ture : COZ 



M and BLW 



25| . as well as one description based on lattice QCD calculations 
26|. The corresponding parameters are presented in Tabled One notices that the param- 
eters r_|_ und r_ in the BLW model and in the lattice calculations are nearly comparable, 
whereas the overall normalization /^r is about 2/3 smaller for the lattice DA as compared 
with the description of the BLW model. In contrast to the BLW model and lattice calcu- 
lations, the parameters r_|_ and r_ are about 3 times larger in the COZ description of the 
nucleon DAs. Below, we will provide calculations using the first two models, COZ and BLW. 
The results following from the lattice calculations can easily be approximated by scaling the 
BLW results. All parameters from the Table [T] have been evolved with leading logarithmic 
accuracy. 



Using the parametrization of Eq. ( 128|1 . the convolution integrals can be computed and 
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yield : 



27r 



2 "J OlemOis 



9 



4?7?7ln ( — ) + (?7 — ?7)(1 — 77777?^) 



+ $2 



(29) 



- ?>{rjrjf{rj - rj) Iv? f ~J + w(l - ^'^VV) [j^) '^^'^ ~ 'r])r]r]{l - r/r/vr^) + |(?7 - ??) 



m 



N 



X 
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2(<l>i + 2$o) (?7-r/)ln 



+ $5 



2?7?7(1 — 67777) In^ 



7777 In^ ( — ) + 1 — r/r/vr^ 



(30) 



12(7/(27/ - 1) In ( - ) + 1 - 12777/ - 27/7/(1 - 67/7/)7r2 



where the notation $j denotes the following combinations of parameters: 

3 1 1 o 3 o 1 



$0 

$1 



1 2 2 
— \ — r_ r r+ — r , H — r+r_ 

4 2 9^2+ + 3^ 

1 1 2 3 5 2 5 

2 9" 2+ 2+ 6 + 



7 2 11 2 4 

18 ~ 2 + 3 ^ 



IV. RESULTS AND DISCUSSION 



We calculate the 27-contribution to the cross section ^27; which is defined by 



dCc.m. — C^0'c.m.,l7 (1 + ^27) ; 



(31) 



as a function of the cm. scattering angle 9, where the cross section da is given by Eq. (|T0|) 
and the cross section in Born approximation dac.m.A'y by Eq. (fT2l) . 

In order to estimate the relative effect of the 27-exchange in the unpolarized cross-section 
we need the input for the time-like FFs in Born approximation, Gm and Ge- We first start 
with a simple description of Gm, given by 

C 



\GM{q')\ 



(32) 
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-1 -0.5 0.5 1 -1 -0.5 0.5 1 

COS cos e 

FIG. 2: Relative 27-contribution to the cross section at (7^ = 6 GeV^ as a function of cos 9. The 
left (right) panels correspond to the calculations using the COZ (BLW) model for the proton DAs 
respectively. The bands describe the contribution for different values of 6Ge given by —6Gm < 
6Ge < SGm- 

with A = 0.3 GeV and with C a free fit parameter. Furthermore, for the first FF parame- 
terization, we use the assumption \Gm\ = \Ge\ and neglect the imaginary part of the FFs 
in our calculation. 

Since the value of SF2 is unknown, we make an estimate of 6Ge using a simple model 

dGE^XSGM, (33) 

where A is a numerical parameter with — 1 < A < 1. 

In our numerical calculations we fix the scale of the running coupling to be /i^ = 0.6 
due to the observation that the scale of the QCD strong coupling is usually lower than the 
value of q^. 

In Fig. [2] we show the relative 27-contribution to the cross section for s = 6 GeV^ as 
a function of cos 9 using the two different models for the proton DAs. The shaded region 
corresponds to the variation of the parameter A in Eq. (!33|) . We observe that for both 
models the relative effect is smaller than 1%. In both cases the angular dependence is 
similar, whereas the COZ model leads to a contribution which is twice as large as when 
using the BLW model. 

The dependence of the 27-contribution on the momentum transfer is shown in Fig. El 
The absolute value of the correction is increasing with g^, but this growth is logarithmic and 
can not change the effect quantitatively. 
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FIG. 3: The magnitude of the relative 27-exchange contribution as a function of for cos = 0.9. 
The left (right) panels indicate the contribution calculated with the COZ (BLW) model respectively. 
The bands describe the contribution for different values of 5Ge given by —6Gm < ^Ge < SGm- 



-0.8 



-1 



— 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

|GeI = IGyl, = Gy* 
impr. log-fit 
IJL fit 




-0.5 





cos9 



0.5 



FIG. 4: Relative 27-exchange contribution for = 6 GeV^ and 5Ge = calculated with the 
COZ model using different parametrizations of the FFs. Solid black curve: using the assumptions 
\Ge\ = \Gm\ and Gm = G\j (purely real FF). Dashed blue curve: improved fit of Ref. [^I involving 
logarithms for F2/F1. Dotted red curve: two-component fit of the FFs as represented in [2]. 



In addition we use two further models to parametrize the FFs entering the 17-amplitude 
and compare the results with the contribution we obtained using the simple fit given by 
Eq. f l32|) . Following we first consider an improved fit of the ratio F2/F1, which includes 
logarithmic corrections to the power law fall-off expected from QCD. Furthermore we use 



a two-component fit of the FFs of Ref. 



28 



The results are demonstrated in Fig. HJ 



where 62^ has been calculated for q"^ = 6 GeV^ with the COZ model for the nucleon DAs. 
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All parametrizations lead to a similar behavior of the 27-contribution with respect to the 
scattering angle and to comparable quantitative results. 

Furthermore we also calculated the single spin asymmetry Ay, which is given by Eq. (IT^ . 
or the two different parametrization of the FFs in Born approximation mentioned above, 
. As discussed in Sec. IIIIl we obtain a small imaginary part of the 27 amplitudes SGm 



27 



28 



and s/mj^F^ in this approach, therefore the 27 contribution to the single spin asymmtry 
mostly results from the interference of the real part of the 27-amplitudes and the imaginary 
part of the FFs Ge and Cm- We find that the relative contribution to Ay is small as well, 
the effect is of order of about 1%. 



V. CONCLUSIONS 

In this work we provided a first calculation of the 27-exchange corrections in pp — )■ e+e~ 
using a pQCD factorization approach. We obtain a small contribution to the cross section of 
< 1%, in the momentum transfer range 5-30 GeV^. The small 27-effect makes it challenging 



to observe such effects in unpolarized cross section measurements, e.g. by PANDA, see |29 |. 
Since the value of the 27-contribution is sensitive to the choice of the DAs, a measurement 
of the process would allow to probe proton and antiproton DAs. 
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